A Double Layer Electromagnetic Cloak And GL EM Modeling 
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In this paper, we propose a novel electromagnetic (EM) cloaking structure that consists of two 
annular layers between three spherical shells and model its performance numerically and theoretically 
by using a Global and Local EM field (GL) method. The two annular layers contain distinct cloaking 
materials: the outer layer provides invisibility; the inner layer is fully absorbing. The cloaking 
materials are weakly degenerative. The wavefield from an EM source located outside the cloak 
propagates as in free space outside the outer shell, never be disturbed by the cloak and does not 
penetrate into the inner absorbing layer and concealment. The field of a source located inside inner 
layer or the cloaked concealment is completely absorbed by the inner layer and never reaches the 
outside of the middle shell. Moreover, the EM wavefield excited in concealment is not disrupted 
by the cloak. There exists no Maxwell EM wavefield can be excited in a single layer cloaked 
concealment which is filled by normal material. Moreover, we find a negative dielectric and positive 
susceptibility metamaterial to fill into the concealment, such that the interior EM wave propagates 
from the concealment to free space through the single layer cloak. Therefore, the double layer cloak 
is important for complete invisibility. Numerical simulations and theoretical analysis verifying these 
properties are performed using the GL EM modeling method that we described in this paper. 

PACS numbers: 13.40.-f, 41.20.-q, 41.20.jb,42.25.Bs 



I. INTRODUCTION 

We propose a novel kind of electromagnetic (EM) dou- 
ble layer cloaking structure (a "Double Layer EM cloak" ) 
that consists of two annular layers between three spher- 
ical shells, Ri < r < R2 and R2 < r < R3. Distinct 
metamaterials are situated in the two layers: the mate- 
rial in the outer layer has properties that provide invisi- 
bility, while material in the inner layer is fully absorbing, 
which can also be useful for making a complete absorp- 
tion boundary condition. We analyze the performance 
of the double layer cloak using a Global and Local EM 
field modeling "GL method" in time domain that we de- 
veloped in this paper, its versions in frequency domain 
have been developed in our other papers [1-3] . Pendry et 
al. [4] proposed a single layer EM cloak material by us- 
ing a coordinate transformation and ray tracing in 2006. 
In their cloak, the ray being bending and re-direction 
around central sphere object and cannot penetrate into 
it. The cloak device like the vacuum and does not dis- 
turb exterior wave field. Later papers studied 2-D plane 
wave propagation through a cloaking layer are published 
in Physical review etc. Journals; Mie scattering analyt- 
ical method for the sphere cloak is studied in [5]; Nu- 
merical methods, including finite-difference time-domain 
method in [6] and finite-element method in [7] are pre- 
sented in which many relative research work papers are 
cited. Most researchers paid attention for exterior EM 
wave field propagation through the cloak. Two authors 
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have studied the cloak's effect on sources located inside 
concealed region [8] and [9]. The authors of [8] stated 
that "when these conditions are over-determined, finite 
energy solutions typically do not exist." They wrote that 
the single layer cloak is insufficient and proposed a dou- 
ble coating to solve this problem. However, there also is 
disputing for this problem. In many numerical simula- 
tions using the GL method, we discovered and verified 
cases where a field satisfying Maxwell's equations can- 
not be excited by a local source inside of the single layer 
cloaked concealment which is filled with normal materi- 
als. For the above case of EM field excited inside the con- 
cealment, our simulations are either divergent or become 
chaotic when the field propagates to the inner boundary 
of the single cloaking layer. These simulations remind 
us to prove that there exists no Maxwell EM wavefield 
can be excited by source inside of the single layer cloaked 
concealment which is filled by normal material. The the- 
oretical proof of the statement is proved in this paper 
by using he GL method. Moreover, we find a negative 
dielectric and positive susceptibility metamaterial to fill 
into the concealment, such that the interior EM wave 
propagates from the concealment to free space through 
the single layer cloak. Therefore, the double layer cloak 
is important for complete invisibility. 

The double layer cloak proposed here overcomes these 
difhculties described above by situating a second, absorb- 
ing layer inside the outer cloaking layer. The wavefield 
excited from an EM source located outside the cloak 
propagates as in free space outside the outer shell and 
never be disturbed by the cloak; it propagates around and 
through the outer layer cloak material and does not pen- 
etrate into the inner layer and innermost concealment. 
The field of a source located inside inner layer or the 
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cloaked concealment is completely absorbed by the inner 
layer and never reaches the outside of the middle shell. 
Un Maxwaell physical behavior at the inner boundary of 
single cloaking layer does not appear in our simulations 
of the double layer cloak. The EM wavefield excited in 
concealment is not disrupted and has no reflection by 
the double layer cloak. In which the EM enviroment in 
the concealment is keeped to be normal and has no be 
changed. As a result, a double layer EM cloak appears 
to be more thoroughly concealing and more robust than 
single layer structures. Moreover, the metamaterials sit- 
uated in the double layer cloak are weakly degenerative. 

The analytical method and numerical method for phys- 
ical simulation have been developed separately in history. 
The GL method consistently combines both analytical 
and numerical approaches. The GL method does not 
need to solve large matrix equations, it only needs to 
solve 3x3 and 6x6 matrix equations. Moreover, the 
GL method does not prescribe any artificial boundary, 
and does not need a PML absorption condition to trun- 
cate the infinite domain. The Finite Element Method 
(FEM) and Finite Difference Method (FDM) have nu- 
merical dispersions which confuse and contaminate the 
physical dispersion in the dispersive metamaterials. The 
frequency limitation is also a difficulty of FEM and FDM. 
The ray tracing is used for wave direction only that is not 
suitable to study full wave field propagation. 

The GL method is a significant scattering process 
which reduces the numerical dispersion and is suitable to 
simulate physical wavefield scattering in the materials, in 
particular, for dispersive materials. Born approximation 
is a conventional method in the quantum mechanics and 
solid state physics. However, it only has one iteration in 
the whole domain which may not be accurate in the high 
frequency and high contrast materials. The GL method 
divides the domain as a set of sub-domains or sub-lattices 
which arc as small as accurate requirement. The Global 
field is updated by the local field from the interaction be- 
tween the global field and local sub-domain materials suc- 
cessively. Once all sub-domain materials arc scattered, 
the GL field solution obtained turn out to be more ac- 
curate than the Born approximation. Moreover, the GL 
method can be mesh-less, including arbitrary geometry 
sub-domains, such as rectangle, cylindrical and spheri- 
cal coordinate mixed coupled together. It is full parallel 
algorithm. The GL method advantages help overcome 
many historical obstacles described in detail in [1]. The 
theoretical foundation of the GL method is described in 
the paper [2]. In particular, for the radial dependent 
cloak materials, the GL method can be reduced to the 
system of the one dimensional GL processes in radial in- 
terval [Ri , R3] by using spherical harmonics expansion. 
We have used the GL modeling [1-2] and inversion tech- 
nique in [3] to simulate many radiaul cloak metamateri- 
als, nanometer materials, periodic photonic crystals etc. 
These simulations show that the GL method is fast and 
accurate. 

The introduction is described in Section 1. The rest 



of this paper is organized as follows: Section 2 describes 
the geometry and material properties of the proposed 
double layer cloak. Section 3 contains a brief descrip- 
tion of the EM integral equations and their solution by 
the GL method. These equations are used in Section 4 
to prove certain theorems about double layer and single 
layer cloaks. Section 5 shows numerical simulations of 
the double layer cloak described in Section 2. The ad- 
vantages of The GL cloak and GL modeling are presented 
in Section 6. Section 7 presents our conclusions. 



II. DOUBLE LAYER ELECTROMAGNETIC 
CLOAK 

In this section, by simulations using GL EM modeling 
and inversion in time domain which is presented in next 
section, we propose a novel electromagnetic (EM) cloak- 
ing structure that consists of two annular layers between 
three spherical shells. An inner layer EM anisotropic 
cloak metamaterial is situated in the inner annular layer; 
Outer layer EM cloak material is situated in outer layer. 



A. Inner Layer EM Anisotropic Cloak Material 

On the inner sphere annular layer domain, fJ/ = 
{r : Ri < r < R2} , we propose an anisotropic material 
as follows, 

[D]j = diag [ei,/Xj] , 

Ei = diag [er,i, £e,ij ^O) 

fii = diag [nr,i, M^.d Mo, ^ ^-^^ 

s,, log 

The flj is called as the inner layer cloak, the material, 
[D]j = diag [si, in (1), is the anisotropic inner layer 
cloak metamaterial. The field of a source located inside 
inner layer or the cloaked concealment is completely ab- 
sorbed by the inner layer and never reaches the outside 
of the middle shell, also, the EM wavefield excited in con- 
cealment is not disrupted by the cloak. The inner layer 
metamaterial, in equation (1), cloaks outer space from 
the local field excited in the inner layer and concealment, 
which can also be useful for making a complete absorp- 
tion boundary condition to truncate infinite domain in 
numerical simulation. 



B. Outer Layer Anisotropic Cloak Material 

Let the outer sphere annular domain Qq = 
{r : R2 < r < R3} he the outer layer cloak with the fol- 
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lowing anisotropic outer layer cloak material, 

[D]o = diag[eo,fi.o] , 

£o = diag [£^,0, ^o.o- £0,0] 

yUo = diag [/U^.o, M0.O. M0.O ] 

^ -„ _fl3i!^\^3 (2) 

se.o = fM).o = s<j>,o = IJ'4>,o 
= . . ^ 

Outer layer cloak provides invisibility, does not disturb 
exterior EM wave field, and cloaks the concealment from 
the exterior EM wavefield. 

C. GL Double Layer Cloak 

The inner cloak 0/ domain and outer cloak CIq domain 
are bordering on the sphere annular surface r = i?2- We 
assemble the as the inner sphere annular domain and 
Slo as the outer sphere annular domain and make them 
coupling on their interface boundary annular surface r = 
i?2 as follows, 

^GL = ni[jno 

= {r:Ri<r<R2}[j{r-R2<r<R3} (3) 
= {r : i?i < r < i?3} , 

and situate the Double Layer anisotropic dielectric and 
susceptibility tensors [D]q]^ on the Ogl as follows. 

The double layer EM cloak materials [D]j- = diag [ei,jli] 
in (1) on the and outer layer cloak material [D]q = 
diag [so, Mo] in (2) on Qq are assembled into the doubled 
layer cloak material [-D]^^ on the domain CIgl- The do- 
main Qgl with the mctamaterial [-D]^^ in (4) is called as 
the double layer cloak which is named GL cloak (GLC). 
GL cloak means that the outer layer has invisibility, never 
disturb exterior wavefield and cloaks the Local conceal- 
ment from the exterior Global field, while the inner layer 
cloaks outer Global space from interior wavefield excited 
in the Local concealment. The properties of the double 
layer material arc studied by GL method simulations and 
analysis in next sections. 



as the outer layer in (4) to construct a GLPS double layer 
cloak. 



E. A Single Negative Refraction Metamaterial 

For studying the EM wavefield excited by source inside 
a concealment which is cloaked by a single layer cloak, 
we propose and situate a following negative dielectric and 
positive susceptibility metamaterial, 

[^]gn = diag [egn,flgn] , 
Egn = diag [e^, ee, e^,] eq, 
p-gn = diag [fir, fie, m] Mo, 

inside the concealment, where r < Ri and r > ro > 0. 
Simulations by the GL EM modeling show that the EM 
wavefield excited by sources inside the c;onceaImcnt with 
metamaterial [Dgn] in (6) will propagate from the con- 
cealment to free space through the single layer cloak. 
Therefore, the double layer cloak is necessary for suffi- 
cient invisiblity cloaking. 



III. 3D GLOBAL AND LOCAL EM MODELING 

In this section, 3D EM integral equations are presented 
briefly. Based on the integral equations, we propose a 3D 
Global and Local EM field modeling. 



A. 3D ELECTROMAGNETIC INTEGRAL 
EQUATION 



D. GLPS Double Layer Cloak 

To compare the properties between the GLO cloak 
and PS cloak, we use the GLI cloak, in equation (1), 
as inner layer and the PS cloak in [4] with the following 
anisotropic dielectric and susceptibility, 

£ps = diag [er,p, £e,p, £<)>,p] £0i 
p.ps = diag [fir^p, iJ'0,p, M0.p] Mo, 



The 3D EM integral equations in frequency domain 
have been proposed in early papers [1] and [2]. In this 
paper, we present the EM integral equations in time 
domain. The integral equations can be derived from 
Maxwell's equations for the electric and magnetic fields, 
E{r,t) and H{r,t): 





' E{r,t)' 






(5) 


H{r,t) 




Hb{r,t) 



n 



Eb{r',t) 
Hb{r',t) 



(7) 
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and 



' E{r,t) ' 






Hir,t) 







E{r',t) 
H(r',t) 



dr'. 



(8) 



In these equations, E and H are the total electric and 
magnetic fields in a given medium with prescribed EM 
sources; Eb and Hb are the fields in a background or ref- 
erence medium with the same sources (the background 
fields are also called "incident" fields); and G 



J,M 
E,H,b 



equation (8) is the EM Green's tensor for the background 
medium, 



Ei^{r',r,t) E^{r',r,t) 



(9) 



The fields E, H, Ei,, and Hi, are all vector quantities. 
The elements of the partitioned Green's tensor are 3x3 
tensors (matrices). For example, a point electric current 
source J in the background medium produces the elec- 
tric field Ejj , while a point magnetic source M produces 
the electric field E^^ . Similar interpretations hold for 
Hjl and H^^ . Also, S [D] is the electromagnetic material 
parameter variation matrix. 



6[D] = 

SDii = (a(r) 
SD22 = ifiir) 



6D11 
6D22 

a,/) + (e-(r)-£,/)f. 



(10) 



9t' 



SDii and SD22 are a 3 x 3 symmetry, inhomogeneous di- 
agonal matrix for the isotropic material, for anisotropic 
material, they are an inhomogeneous diagonal or full ma- 
trix, / is a 3 X 3 unit matrix; a{r) is the conductivity 
tensor, e(r) is the dielectric tensor, p,{r) is susceptibility 
tensor of the full medium, which can be dispersive (i.e., 
its properties can depend on the angular frequency w); df, 
is the conductivity, Sb is the permittivity, fXb is the perme- 
ability in the background medium; f2 is the finite domain 
in which the parameter variation matrix S [D] ^ 0; the 
(e(r) — £bI)E is the induced electric polarization; and 
(/x(r) — HbI)El is the induced magnetization. Finally, *t 
is convolution with respect to time, t. 



(3.2) When fc = 0, let Eo{r,t) and HQ{r,t) are the an- 
alytical global field, E;^{r', r, t), {r' , r, t), E^ {r' , r, t), 
and {r' ,r,t) are the analytical global Green's ten- 
sor in the background medium. By induction, suppose 
that Ek-iir,t), Hk-iir,t), E^^^ir' ,r,t), H(_^{r' ,r,t), 
E^'_T^{r' ,r,t)^ and H^'^-^^{r' ,r,t) are calculated in the 
(fc — 1)*'* step in the subdomain flk-i- 

(3.3) In {i^fc}, upon substituting Ek-i{r,t), 
Hk-i{r,t), Ei_,{r',r,t), Hi^,{r' ,r,t), E^L,{r' ,r,t), 
and H^'^-^{r' ,r,t) into the integral equation (8), the EM 
Green's tensor integral equation (8) in flk is reduced 
into 6x6 matrix equations. By solving the 6x6 matrix 
equations. 



we obtain the Green's tensor field Ej:{r' , r, t), 



H^ir',r,t),Ei'{r\r,t). 



and H^^{r',r, t). 



(3.4) According to the integral equation (7), the elec- 
tromagnetic field Ek (r, t) and {r, t) are updated by 
the interactive scattering field between the Green's ten- 
sor and local polarization and magnetization in the sub- 
domain 0/j as follows. 



" -Bfe(r,t) ■ 




' Ek-i{r,t) ' 






Hk-i{r,t) 



+ 1 



tS[D{r') 



Ei{r\r,t) Hi{r',r,t) 
E^^{r',r,t) H^\r\r,t) 

Ek-i{r',t) 
Hk-i{r',t) 



(11) 



dr' 



(3.5) The steps (3.2) and (3.4) form a finite iteration, 
k — 1, 2, • • • ,N, the En (r, t) and Hn (r, t) are the elec- 
tromagnetic field of the GL modeling method. The GL 
electromagnetic field modeling in the time space domain 
is short named as GLT method. 

The GL EM modeling in the space frequency domain 
is proposed in the paper [2[ , we call the GL modeling in 
frequency domain as GLF method. 



IV. INTERACTION BETWEEN THE EM 
WAVE FIELD WITH CLOAKS 

A. Theoretical Analysis Of Interaction Between 
The EM Wave Field With The GL Double Layer 
Cloaks 



B. 3D GL EM MODELING 

Although the two integral equations (7) and (8) are 
equivalent, there are numerical advantages in using both 
sets simultaneously. Based on the integral equation, the 
GL EM modeling in the time domain is proposed in this 
section. The 3D GL EM modeling and integral equations 
in the frequency domain have been developed in early 
papers [1-2]. 

(3.1) The domain is divided into a set of N sub 

N 

domains, {rifc}, such that = IJ flk- The division can 

fe=i 

be mesh or meshless. 



Theoretical analysis of the interaction between the EM 
wave with GL cloaks are presented in this section. 

Statement 1: Let domain ^Igl, equation (3), and 
anisotropic metamaterial Dgl, equation (4), be a "GL" 
double layer cloak, and let the concealed central sphere 
|r| < i?i and the region outside the cloak |r| > i?3 be 
filled with a normal electromagnetic material with per- 
mittivity and permeability, e—eh,fJ- — fJ-b- We claim the 
following: (1) provide EM wavefield is excited by sources 
located inside the concealment of GL cloak, |fs| < 
the EM wavefield never be disturbed by the cloak; (2) 
provide EM wavefield is excited by sources located inside 
concealment or the inner layer of the cloak, \fs\ < R2, 
the EM wavefield is absorbed in the inner layer Qj and 
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docs not extend outside the boundary r = i?2- (3) pro- 
vide the wavefield is excited by sources located outside 
of the cloak, \ fs\ > R3, the EM wavefield propagates as 
if in free space, undisturbed by the cloak; (4) provide 
the EM wavefield is excited by sources outside cloak or 
inside outer layer, \rs\ > R2, the EM wavefield does not 
penetrate into the inner layer and the concealed region, 
1^ < R2. 



B. No Maxwell EM Wavefield Can Be Excited By 
Nonzero Local Sources Inside The Single Layer 
Cloaked Concealment Fielled With Normal 
Materials 

Statement 2: Suppose that a 3D anisotropic inho- 
mogeneous single layer cloak domain separates the whole 
space into three sub domains: one is the single layer cloak 
domain fldk with the cloaking material; another is the 
cloaked concealment domain flconi with normal EM ma- 
terials; and the third is free space outside the cloak. If 
the Maxwell EM wavefield excited by a point source or 
local sources outside the concealment Slconi vanishes in- 
side the concealment ilconi , then there is no Maxwell EM 
wave field excited by the local sources inside the cloaked 
concealment flconi- 

The Maxwell EM wavefield is an EM wave field that 
satisfies the Maxwell equation and appropriate interface 
boundary conditions. We use an inverse process to prove 
the statement 2 as follows: Suppose that there exists a 
Maxwell EM wavefield excited by the local sources in- 
side the concealment, filled with the normal materials, 
the wavefield will satisfy the Maxwell equation in all of 
the space, including in the anisotropic inhomogeneous 
cloaking layer flcik and in the concealment, flconi] and 
also satisfies boundary conditions on the interfaces, Si 
and 52. Si is the interface boundary between the cloak- 
ing region Qcik and the concealment domain ilconi, and 
also is the inner boundary surface of the cloaking domain 
^cik- S2 is the interface boundary between the cloaking 
region ilcik and free space outside; it also is the outer 
boundary surface of the cloaking region fldk- 

Let Rc = - flcik U ^coni, Rd = R^ - ^coni, and by 
integral equation (7), the EM wave field satisfies 



becomes to 



' E{r,t)' 




■ Eb{r,t) ' 


_H{r,t) 




_H,ir,t)_ 



+ 



I Gi%{r',r,t)nS[D] 



Eb{r',t) 
Hk{r',t) 



dr', 



(12) 



where G'^^^ {r' , r, t) is the Green's tensor, its components 
E-', H-'. E^, and H^^{r',r,t) are the Green's func- 
tion on flcik[j^coni[J Rc, excited by the point impulse 
sources outside the concealment, r ^ Rd- By the assump- 
tions, (r', r, t) exists on ^cik U ^coni U and when 
r' € flconi, G"^^ (r', r, t) = 0. The integral equation (12) 



' E{r,t) ' 




■ Eb{r,t) ' 


H{r,t) 




Hb{r,t) 



+ I G-^/;^{r',r,t)HS[D] 



Eb{r\t) 
Hb{r',t) 



dr'. 



(13) 



Consider the Maxwell equation in i?^, the virtual 
source is located r, r £ Rj, and the point source is lo- 
cated rs, the variable in the following equations is r', 
rs G ^coni and Vg ^ Rd, we have 



Vx 



-Vx 



Gi%{r',r,t) 



(14) 



{r',r,t)+I6{r',r)S{t), 



and 



-Vx 
= [Db] 



Vx 

Eb' 
Hb 



Eb 
Hb 

{r',rs,t), 



{r',rs,t) 



(15) 



Convolving equation (14) with [Eb {r',t) , Hb (r'.t)], and 
equation (15) with (r', r, t), subtracting the second 

result from the first, integrating by parts, and making 
some manipulations gives 



Ebir,t) 
Hb{r,t) 

+ I G'^%ir',rA)nS[D 



= §G'i'%{r',r,t)®t 

Si 



Ebir',t) 
Hb{r',t) 

Eb{r',t) 
Hb{r',t) 



dr' 



(16) 



dS, 



where (8)t is cross convolution. From the assumption 
of the statement 2 that "the Maxwell EM wavefield ex- 
cited by a point source or local sources outside the con- 
cealment ilconi vanishes inside the concealment riconj" j 
and virtual source r is located outside the concealment, 
r G Rd, if r' G ftconi, G'^'^ {r',r,t) = 0. By continuity, 
when r' G Si, wc have G'^^fj {r' ,r,t) = 0. Continuity of 
G^E^ (r',r, f) at interface implies that the term in right 
hand side of (16) vanishes, giving 



+ 



Eb{r,t) 
Hb{r,t) 



+ / G^;^,(,(r',r,i)*,5[D] 



E{r',t) 
H{r',t) 



dr' = 0. 



(17) 



Substituting equation (17) into the integral equation 

(13) gives 



Eir,t) 
H{r,t) 



0. 



(18) 
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From continuity condition on the EM wave field, we ob- 
tain the fohowing over vanishing condition on the bound- 
ary of the concealment Slconi, 



E{r,rs,t) 
H{r,rs,t) 



0, 



(19) 



where rs denotes point source location inside of the con- 
cealment, r is EM field receiver point, r G Si. Because 
the EM wavefield is excited by local sources inside of 
the concealment domain ^coni, it satisfies the following 
Maxwell equation with variable r', 

Vx 



-Vx 

= [Dconl] 





■ E ' 




H 



E 
H 



ir',r„t) 
{r',r„t) + Q{r',r„t). 



(20) 



where [Dconi] = diag[er6Ql, j-iri-J-al] (d/dt) with the 
normal EM material parameters, > 1 and /ir > 1 are 
relative EM parameters, Eq is basic permittivity and /io is 
basic permeability, / is 3 x 3 unit matrix, rg S i^coni is the 
local source location, Q{r' , rs,t) is the nozero local source 



inside the concealment, O.conl- 
Greens tensor which satisfies 



Let G 



J,M 

E,H,conl 



Vx 
-Vx 

+IS{r',r)6it) 



{r',r,t) 



{r',r,t) be 



(21) 



Convolving [E (r' , t) , H {r' , t)] with (21), and 
G'^'% ^g^i{r' ,r,t) with (20), subtracting the second 
result from the first, taking the integral to flconl, and 
use integration by parts and other manipulations, we 

have 



E{r,rs,t) 
H{r,rs,t) 

■ I Gi%coniir',r,t)HQ{r',r,,t)dr' 



(22) 



§ G 



E,H,conl 



(r',r,t) (g)t 



E{r',rs,t} 



dr' 



(81* denotes the cross convolution, and dilconi = Si. Be- 
cause of the over vanishing condition (19), 



E{r, r,,f) 
H (r, rs,t) 



= 0, 



Si 



we have 



E{r,rs,t) 
H{r,rs,t) 



i{r',r,t) *tQ{r',rs,t)dr'. 



(23) 



Because G'^^^^^^^i{r' ,r,t) ^ and Q{r',rs,t) ^ 0, so, 



From the continuity of the EM wavefield, the nonzero 
EM wave field (24) results that 



Eir,rs,t) 



7^0. 



(25) 



The EM wavefield is nonzero on the boundasry in 
(25) that is an obvious contradiction with the same EM 
wavefield is zero on the boundasry Si in (19). Therefore, 
we proved statement 2. Result (24) can also be derived 
more simply and obvious from the integral expression 
(23). Let the source be a point impulse current source 
with polarization direction x, i.e.. 



(23), gives 



Q{r,rs,t) = 


--S{r 


(26) and Sr 


= ] 


' E{r,r,,t) ■ 




_H{r,rs,t) 





Ei{r,r,,t) 



Exx {r. rs , t) 
E'xy {i\ rg , t) 
Exz{r,rs,t) 



(26) 

1.0 into the 
(27) 

(28) 



Exxir. rs.t) 

1 3^ S(t-,/sJI\r-r,\) 

Stz^c dx^ \r—rs\ 
1 02 S[t-^\r-rs\) 



+ 8^2/^9*2 \r-rs\ 

It is obvious that when r G Si 

Exxir,rs,t)\^^g^ ^ 0. 



(29) 



(30) 



The electric intensity field Exx{r,rs,t)\j.^g^ ^ in (30) 



and Exxir,rs,t)\ 



0. in (19) are an obvious contra- 



diction. The Statement 2 is true. 



V. SIMULATIONS OF THE EM WAVE FIELD 
THROUGH THE DOUBLE LAYER CLOAK 

A. The Model of The Double Layer Cloak 

The full 3D simulation model is a unit cube, [0.5, 0.5]'^, 
centered on the origin, discretized on a 201^ mesh, with 
uniform mesh spacing of 0.005m. The EM wavefield is 
excited by a point electric source: 



s{r,rs,t) =6{r- rs)d{t)e, 



(31) 



at location where e is the (unit) polarization vector, 
the time step dt = 0.3333 x 10~^° ; the frequency band 
is from Q.Q^GHz to 15GHz. The shortest wavelength is 
about 0.02m. The GL double layer EM cloak consists 
of inner and outer annular regions, floL = ^i[J^o, 
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equation (3), with centers at the origin, which is situ- 
ated by proposed anisotropic metamaterial -Dgl, equa- 
tion (4), the concealed central sphere |r| < i?i and the 
region outside the cloak |r| > i?3 is filled with a nor- 
mal electromagnetic material with basic permittivity and 
permeability, s—eh,fJ- = fJ-t- The inner boundary of the 
cloak is Ri = 0.2m; the middle shell boundary between 
the two layers is i?2 = 0.3m; and the outer boundary 
is i?3 = 0.45m. In the simulations, the sphere region 
r < R3 is actually modeled in spherical coordinates, 
{r,6,ip), where 9 is polar angle. The sphere is divided 
180'^ cells. The spherical coordinate grid is superimposed 
on the rectangular grid used to mesh the domain outside 
shell r = R3. Because the cloaking materials are radial 
dependent, the EM modeling is reduced to the system 
of the one dimensional GL modeling by using the sphere 
harmonic expansion. 




jv s hi I I 1 1 I 1 1 I I 1 1 I I 1 1 I I 1 1 I I 1 1 I 1 1 1 I 1 1 1 I 1 1 I I 1 1 I I 1 1 1 1 1 I i i 
41.5 -CA -C.S -0.2 ^.r C- G.1 D.: D.l C.'S n.5 
Horizontal X IMeterS. Made In GLGEO 



FIG. 1: (color online) At time step 18dt, the wave front of the 
First electric wave, Exx,i, propagates inside the concealment 
r < Til. The front of Second EM wave, Exx,2, is located in 
free space, the right and top corner outside of the whole GL 
double layer cloak. 



ESectriG Wavie Era Propjgalkin in Double Layered EM GL-Ckjak 




FIG. 2: (color online) At this moment of the time step 
30dt, the front of the First electric wave, Exx.i, propa- 
gates enter to the inner layer, Ri < r < R2; the front 
of Second electric wave, Exx,2, reaches the outer boundary 
r = f?3 of the GL double layer cloak. 

ElectTK] Wave Exx FrDp|aga1ion rn Dibble Layered EM GL-Cloak 
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■-0.5 -0.4 -C.3 -a.: -D.1 0.1 0.2 0.3 OA C.5 
Hofizomal X (Meter), Made In GLGEO 



B. Simulation I 

In this subsection, the simulation I is presented that 
an inner point source in the concealment and other outer 
source in the free space are used to excite the EM wave 
propagation through the double layer cloak. The point 
impulse current source with polarization direction x, i.e. 

S{r,rs,t)^Sir~rs)Sit)x. (32) 

The first point current source is located inside the 
concealment at (0.0012m, 0.0m, 0.0m), by which the ex- 
cited EM wave is named as First EM wave, its com- 
ponent Exx,i is labeled First electric wave. The sec- 
ond current point source is located in free space at 



FIG. 3: (color online) At the time step 58cit, the 
First electric wave, Exx.i, is propagating inside the inner 
layer, R\ < r < R2, and becomes very slow. The part of 
the front of Second electric wave, Exx,2, has been inside of 
outer layer, R2 < r < R3, and being backward bending 



(0.518m, 0.518m, 0.0) where is in the right and top cor- 
ner outside the double layer cloak. The EM wave by the 
second source is named as Second EM wave. Its com- 
ponent Exx,2, is labeled Second electric wave. Figure 1- 
6 show a series of snapshots of an EM field propagat- 
ing in and around the double layer cloak. The Figure 
1 shows that at time step 18dt, the wave front of the 
First electric wave, Exx.i, propagates inside the central 
sphere concealment r < Ri and never be disturbed by 
the inner layer; The front of Second EM wave, Exx,2, 
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FIG. 4: (color online) At time step 75cit, the 

Second electric wave, Exx,2, is propagating inside of the outer 
layer, R2 < r < R3, and around the shell r — R2 and 
never penetrate into inner layer and concealment, r < R2. 
It does split into the two phases around the shell r = i?2- 
The First electric wave, Exx.i, is propagating inside the in- 
ner layer, Ri < r < R2. 




-o.s -0.4 -0.3 -a.: -D,i D 0.1 0.2 0.3 sa 
Horizontal X {Meter), MatJe In GLGEO 



FIG. 6: (color online) At time step 128dt, the 
Second electric wave, Exx,2 has propagated outside the GL 
double layer cloak, a small part of its wave front is located 
in the left and low corner of the plot frame, which never be 
disturbed. The First electric wave, Exx,i, is still propagating 
inside the inner layer, Ri < r < 7?2. 



ElMfric Ware Ejot Propagafon in Double Layered EM GL-C leak 
In XY Plane UEing 3D GL Modeling in G LGEO 
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FIG. 5: (color online) At time step 9Sdt, one part of front 
the Second electric wave, Exx,2, is propagating inside of the 
outer layer, R2 < r < R3, and around to other side far source. 
It is a little forward bending. The First electric wave, Exx,i, 
is still propagating inside the inner layer, Ri < r < R2. 



is located in free space, the right and top corner out- 
side the double layer cloak. At the moment 30dt in Fig- 
ure 2, the front of the First electric wave, Exx.i, prop- 
agates enter to the inner layer, Ri < r < R2; The 
front of Second electric wave, Exx,2, reaches the outer 
shell boundary r = i?3 of the cloak and never be dis- 
turbed. At the time step 58dt, the First electric wave, 
Exx.ii is propagating inside the inner layer, Ri < r < 




■-D.5 -■;.4 -C.3 -a2. ^.T Q G.I D.2 D.3 v.i D.5 
HoFtzantal X (Metef), Made In GLGEO 



FIG. 7: (color online) In the double layer cloak, inner layer 
is situated with GLI cloak material in (1), outer layer is 
filled with PS cloak material. At the time step 58dt, the 
First electric wave, Exx,i, is propagating inside inner layer; 
The Second electric wave propagates inside outer layer with 
more backward bending and more slow than EM wavefield in 
GLO material (2) in Figure 5. 



R2, and becomes very slow; The part of the front 
of Second electric wave, Exx,2i propagates inside outer 
layer, i?2 < < ^3, and being backward bending. 
The EM wave propagation image snapshot is presented 
in Figure 3. In the Figure 4, at time step l^dt, the 
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FIG. 8: (color online) In the double layer cloak 
GLPS arranged in Figure 7, at the time step 75dt, the 
First electric wave, Exx.i, is propagating inside inner layer; 
The Second electric wave propagates around the shell r = R2 
with more forward bending than EM wavefield in GLO ma- 
terial (2) in Figure 5. 



Elecjtiic Wave Pn^aga^on in Dovble Lawred EM -CkHk 
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FIG. 9: (color online) In the double layer cloak GLPS in 
Figure 7, at the time step 98dt, the First electnc wave, Exx,i, 
is propagating inside inner layer; The Second electric wave 
propagates inside outer layer around to side far source with 
more forward bending and more fast than EM wavefield does 
in GLO material (2) in Figure 5. 



First electric 



Exx,i, is still propagating inside the 



inner layer, Ri < r < R2; The Second electric wave, 
Exx,2, is propagating inside the outer layer cloak, R2 < 
r < R3, and around the shell r = R2 and never penetrate 
into the inner layer and concealment, r < i?2- It does 
split into the two phases around the shell r = R2, the 



front phase speed exceeds the light speed; the back phase 
is slower than the light speed. In the figure 5, at time step 
98dt, one part of front the Second electric wave, Exx,2, 
is propagating inside the outer layer , R2 < r < R3. 
has a little forward bending and never penetrate into 
the inner layer and the concealment, i.e. r < R2. The 
First electric wave, Exx,i, is still propagating inside of 
the inner layer, i?i < r < i?2. At time step 128dt, the 
Second electric wave, Exx,2 has propagated outside dou- 
ble layer cloak, a small part of its wave front is located in 
the left and low corner of the plot frame which is shown in 
the Figure 6, most part of front of the Exx,2 electric wave 
field has been out of the plot frame The exterior EM wave 
outside the cloak never been disturbed by the cloak and 
never penetrate enter its concealment and inner layer. 
At same time step, the First electric wave, Exx,i, is still 
propagating inside the the inner layer, i?i < r < i?2- 
It can be very closed to the interface boundary shell 
r — R2, However, it can not be reached to the boundary 
r = i?2 for any long time. That means that the interior 
EM wavefield excited by source inside the concealment 
is complete absorbed by the inner layer. The inner layer 
metamaterial, in equation (1), cloaks outer space from 
the local field excited in the inner layer and concealment, 
which can also be useful for making a complete absorp- 
tion boundary condition to truncate infinite domain in 
numerical simulation. 
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FIG. 10: (color online) At time step 30dt, one part of the front 
of the Second electnc wavefield, Exx,2, propagates in free 
space with disturbance, its other part of the front has been 
propagating around the shell r — R2; It does not propagate 
into the inner layer and concealment. The First electnc wave, 
Exx.i, is very slow and just starts to propagate in the inner 
layer. 
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FIG. 11: (color online) At time step 68dt, front of the 
First electric wave, Exx,i, is very slow propagating inside 
inner layer R\ < r < R2; The most part of front of 
Second EM wave, Exx,2, propagates in free space with dis- 
turbance, its other small part of front propagates in the outer 
layer around to side far source, but it does not propagate into 
the inner layer. 




-0.4 -[i.i -:.2 -d.i c oi z-.i a 04 d.6 
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FIG. 13: (color online) The PS single layer cloaked con- 
cealment is filled with a permittivity single negative metama- 
terial, The electric wavefield Exx excited by a point source 
inside the concealment propagates out to free space through 
the cloak. 



Elec^oWave Exji Pmpagiiaon in DDubJe Layered EM GL-Cioak 




FIG. 12: (color online) At moment 98dt, one part of the 
front of the First electric wave, Exx,i, propagates enter to the 
concealment; other part of the front is still propagating inside 
of inner layer, i?i < r < 7?2., it does not propagate outside 
shell r — R2; The front of Second electric wave, Exx,2, has 
propagated outside all cloak, only very small part of the front 
is located the left lower corner of the figure frame. 



C. Simulation II 

To compare the properties between the GLO cloak 
and PS cloak, we use PS cloak [4] as outer layer and 
GLI cloak, in equation (1), as inner layer to construct 
a GLPS double layer cloak. The simulations of the EM 



wavefield propagation through the GLPS cloak is pre- 
sented in simulation II. The cloak geometry and mesh 
configuration is described in subsection A. Because the 
GLI cloak, in equation (1), is used as inner cloak of GL 
double layer cloak and GLPS double layer cloak, the 
First electric wave propagation in timestep 30dt, 58dt, 
and 98dt are as same as in Figure 3 and Figure 7; Figure 
4 and Figure 8, Figure 5 and Figure 9, respectively. At 
timesep 58dt, the Second electric wave, Exx,2, propagat- 
ing in PS cloak is more slow and more backward bending 
than the field does in GLO outer layer, in Figure 3. At 
timestep 75dt, the front of Second electric wave, Exx.2i 
propagating in outer layer of GL cloak in Figure 4, is 
better than the field propagation in PS cloak, in Figure 
8. In particular, at timestep 98dt, Second electric wave, 
Exx,2, propagating in layer PS cloak has more foeward 
bending in Figure 9, but the field has a little forward 
bending in GLO outer layer in figure 4. Summary, GLO 
cloak material degeneration is weaker and better than 
the PS cloak. 



D. Simulation III 

The EM wavefield excited by a point source inside the 
inner layer and other point source inside the outer layer 
of GL cloak propagates through the double layer cloak 
that is presented in this subsection. By the first point 
current source at (— 0.18m, —0.18to, 0.0), the excited 
EM wavefield Exx.i is labeled First electric wave. The 
second current point source at (0.258m, 0.258rn, 0.0) 
excites the EM wavefield, Exx.2, which is labeled 
Second electric wave. The Figure 10 shows that 
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at time step 30dt, one part of the front of the 
Second electric wavefield, Exx.2, propagates in free 
spaee with disturbance in [15], its other part of the front 
has been propagating around the shell r = R2; It does 
not propagate into the inner layer and concealment. The 
First electric wcLvej Exx,i^ is very slow and just starts 
to propagate in the inner layer. In the Figure 11, at 
time step 68dt, front of the First electric wave, Exx,i, is 
very slow propagating inside inner layer i?i < r < i?2; 
The most part of front of Second EM wave, Exx,2i 
propagates in free space with disturbance, its other 
part of front propagates inside the outer layer around 
to side far source, but it does not propagate into the 
inner layer. The wavefield propagating at moment 
%Sdt is presented in figure 12. One part of the front of 
the First electric wave, Exx.i, propagates enter to the 
concealment; other part of the front is still propagating 
inside of inner layer, Ri < r < R2-, it does not propagate 
outside shell r = R2; The front of Second electric wave, 
Exx.2, has propagated outside all cloak, only very small 
part of the front is located the left lower corner of the 
figure frame 



E. Simulation IV 



For studying the EM wavefield excited by source in- 
side a concealment which is cloaked by a single layer 
cloak, we propose and situate a novel negative dielec- 
tric and positive susceptibility metamaterial [D]aN, in 
equation (6), inside the concealment, where r < Ri and 
r > ro > 0. The concealment is cloaked by the PS cloak 
[4], in equation (5). The simulation domain and mesh 
configration is described in subsection A. A time har- 
manic point source with polarization in y direction and 
frequency 0.57 x 10^°Hz, 



s(r,rs,t) = 5{r-rs)e'^'y, (33) 



is located in the point r^, (0.0,0.26,0.0). Simulations 
by the GL EM modeling show that, in Figure 13, the 
EM wavefield Eyy excited by the source in (33) inside 
the concealment filled with the metamaterial is propa- 
gating from the concealment to free space through the 
PS single layer cloak. Therefore, the double layer cloak 
is necessary for sufficient invisiblity cloaking. This single 
negative refraction metamaterial, in equation (6), will be 
investigated and detailed presented in next paper. 



VI. ADVANTAGES 

A. The EM GL Double Layer Cloak Is Robust and 
Sufficient For Invisibility 

The figure 10-12 clearly show that at moment of the 
30dt, 68dt, and 98dt time step, the wave front of the 
Second electric wave has propagated outside the cloak 
and go to free space with disturbance. The result re- 
minders us that if only single outer layer cloak Oq, or 
PS cloak is adopted, and there is a little crack loss on 
the inner side of the boundary surface d^o-, some EM 
or current source inside the fio will excite the EM wave 
propagation go out to free space and expose the cloak im- 
mediately. Moreover, the Figure 13 shows that the EM 
wavefield Eyy excited by the source in (33) inside the 
concealment filled with the metamaterial is propagating 
from the concealment to free space through the PS sin- 
gle layer cloak. Such that the single layer cloak complete 
lose the cloaking function. The GL double layer cloak 
overcomes the weakness that is also shown in the figure 
1-12. The wave front of the First electric wave, which is 
excited by a point source inside the inner layer cloak f2/ , 
or inside the concealment, is always propagating inside 
of inner layer cloak f7/ or concealment VLconi and never 
propagate outside of the interface shell r = i?2 • Any EM 
field inside the inner layer or inside the concealment can 
not be propagated outside the shell r = i?2- Therefore, 
the EM GL double layer cloak is robust and sufficient for 
invisibility. 

Using the GL method theoretical analysis, the state- 
ment 2 in section 4 is rigorously proved. It states that 
"there exists no Maxwell electromagnetic wavefield can 
be excited by nonzero local sources inside of the single 
layer cloaked concealment with the normal EM materi- 
als". The invisibility of the single layer cloak and ex- 
istence of Maxwell EM wave field excited by the local 
sources inside its concealment is inconsistent. Provide 
only single outer layer cloak is adopted, the EM field 
excited by local sources inside of its concealment with 
normal materials does not satisfy the Maxwell equation. 
The EM chaos phenomena, which is divorced from the 
Maxwell equation governing, may damage devices and 
human inside the concealment, or may degrade the invis- 
ibility of the cloak. The single layer cloak is not complete 
and unsafe. The Figure 13 shows when a special meta- 
material fills into the concealment, the single layer cloak 
complete lose cloaking function. The GL double layer 
cloak in this paper or double coating in paper [8] are 
necessary for the complete invisibility function. The GL 
double layer cloak is different from the double coating 
in paper [8]. The simulations in Figures 1-12 show that 
the inner layer metamaterial of GL double layer, in equa- 
tion (1), does not disturb and reflect the EM field excited 
by sources inside the concealment and does not change 
the EM environment in the concealment. Therefore, GL 
cloak is safe metamaterial. 
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B. Preqency Band And Geometry 

Many simulations and theoretical analysis by the GL 

method show that the idea EM GL double layer cloak is 
of the invisibility function for all frequencies. However, 
the practical material has some loss. The frequency band 
will be depended on the rate of the material loss. 

When r is decreasing and going to i?2 in outer 
layer, the EM wavefield speed in PS cloak, In equa- 
tion (5), is increasing to infinite with the strong diver- 
gent rate l/(r — R2), but the EM wavefield speed in 
outer layer of GL cloak, in equation (1), is increasing 
with the weak divergent rate l/\/r — i?2- In the side 
near the source, the simulation in Figure 7 shows that 
Second electric wave, propagation is more slow 

and more backward bending inside the PS cloak than the 
field does in outer layer of GL cloak in Figure 3. In the 
side far the source, the simulation in Figure 9 shows that 
Second, electric wawe, 2, propagation is more fast and 
more forward bending inside the PS cloak than the field 
does in outer layer of GL cloak in Figure 5. 

The comparisons between the EM wave propagation 
inside the GL cloak and PS cloak, in Figure 3-5 and Fig- 
ure 7-9, show that the GL cloak has better properties 
to reduce the dispersion and degeneration. Therefore, a 
reasonable frequency band of the GL double layer cloak 
for low loss rate may be obtained in practical fabrica- 
tion. The EM GL double layer cloak can be extended to 
have double ellipsoid annular layer and other double strip 
geometry or use non edclidean geometry [10] in outer 
layer cloak and GLI in (1) in inner layer to make wide 
frequency band double layer cloak in next paper. The 
experiments in this field, the pioneering demonstration 
by D. Schurig et al in [11] and the recent paper by R. 
Liu et al. in [12] may be important help for improving 
cloak model and simulation. As early in 2001, a double 
layer cloth phenomenon has been observed in Lawrence 
Berkely National Laboratory which is published in SEG 
Expanded Abstracts in 2002 [13]. For studying and sim- 
ulating the strange phenomenon and metamaterials, we 
developed the effective GL modeling and inversion. The 
first paper of the GL method and a relative EM "mirage" 
phenomenon have been presented in PIERS 2005 and 
published in the proceeding of PIERS 2005 in Hangzhou 
in [14] and [15]. 

C. Advantages Of The GL Method 

The GL EM modeling is fully different from FEM and 
FD and Born approximation methods and overcome their 
difHculties. There is no big matrix equation to solve in 
GL method. Moreover, it does not need artificial bound- 
ary and absorption condition to truncate the infinite do- 



main. Born Approximation is a conventional method in 
the quantum mechanics and solid physics, however, it is 
one iteration only in whole domain which is not accu- 
rate for high frequency and for high contrast materials. 
The GL method divides the domain as a set of small 
sub domains or sub lattices. The Global field is updated 
by the local field from the interaction between the global 
field and local subdomain materials successively. Once all 
subdomain materials are scattered, the GL field solution 
is obtained which is much more accurate than the Born 
approximation. GL method is suitable for all frequency 
and high contrast materials. 

Moreover, the GL method can be meshless, including 
arbitrary geometry subdomains, such as rectangle, cylin- 
drical and spherical coordinate mixed coupled together. 
It is full parallel algorithm. The GL EM method consis- 
tent combines the analytical and numerical approaches 
together and reduced the numerical dispersion and nu- 
merical frequency limitation. The GL method has dou- 
ble capabilities of the theoretical analysis and numerical 
simulations that has been shown in this paper. Because 
the cloak metamaterial in (4), (5), and (6) are radial 
dependent, the reduced one dimension GL EM model- 
ing system by using the sphere harmonic expansion are 
developed for the cloak simulations. The 3D GL simula- 
tions of the EM wave field through the single and mul- 
tiple sphere, cylinder, ellipsoid, and arbitrary geometry 
cloaks in single layer and double layer cloaks show that 
the GLT and GLF EM modeling are accurate, stable and 
fast. The 3D and 2D GL parallel software are developed 
and patented by GLGEO. 



VII. CONCLUSIONS 

Many simulations by the GL modeling and theoretical 
analysis verify that the EM GL doubled cloak is robust 

and safe cloak and has complete and sufficient invisibility 
functions. Its concealment is the normal electromagnetic 
environment. The outer layer of the GL double layer 
cloak has the invisible function, its inner layer cloak has 
fully absorption function. The GL method is an effec- 
tive physical simulation method. It has double c;apabil- 
ity of the theoretical analysis and numerical simulations 
to study the cloak metamaterials and wide material and 
Field scattering in physical sciences. 
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